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ABSTRACT 

We have used Kepler photometry to characterize variability in four radio-loud active 
galactic nuclei (three quasars and one object tentatively identified as a Seyfert 1.5 
galaxy) on timescales from minutes to months, comparable to the light crossing time 
of the accretion disk around the central supermassive black hole or the base of the 
relativistic jet. Kepler's almost continuous observations provide much better temporal 
coverage than is possible from ground-based observations. We report the first such 
data analyzed for quasars. We have constructed power spectral densities using 8 Kepler 
quarters of long-cadence (30-minute) data for three AGN, 6 quarters for one AGN and 
2 quarters of short-cadence (1-minute) data for all four AGN. On timescales longer than 
about 0.2-0.6 day, we find red noise with mean power-law slopes ranging from -1.8 to 
-1.2, consistent with the variability originating in turbulence either behind a shock or 
within an accretion disk. Each AGN has a range of red noise slopes which vary slightly 
by month and quarter of observation. No quasi-periodic oscillations of astrophysical 
origin were detected. We detected several days-long flares when brightness increased 
by 3% - 7% in two objects. No flares on timescales of minutes to hours were detected. 
Our observations imply that the duty cycle for enhanced activity in these radio-loud 
AGN is small. These well-sampled AGN light curves provide an impetus to develop 
more detailed models of turbulence in jets and instabilities in accretion disks. 

Subject headings: accretion disks - black hole physics - galaxies: active - galaxies - 
quasars: general - galaxies: Seyfert 

1. Introduction 

One of the best ways to probe the extremely small regions from which the bulk of the energy 
in AGN is emitted is through the study of their variability in different bands. The immense 

^Space Science Institute, 4750 Walnut Street, Suite 205, Boulder, CO 80301 

^Department of Physics, The College of New Jersey, PO Box 7718, Ewing, NJ 08628, USA 

^Jet Propulsion Laboratory, Mail Stop 321-100, 4800 Oak Grove Drive, Pasadena, CA 91109 



- 2 - 



13 powers, non-thermal spectra, and rapid variability detected across the electromagnetic spectrum 

14 that characterize the class of blazars, i.e., Flat Spectrum Radio Quasars (FSRQs) and BL Lac 

15 objects, can only be understood within the framework of matter flowing inwards through accretion 

16 disks (ADs) onto supermassive black holes (SMBHs). The bulk of this emission emerges from 

17 within several gravitational radii of the SMBH; part of the energy often is channelled outward via 
IB relativistic jets, producing radio- loud AGN. That emission is magnified by Doppler boosting when 
19 the jet is pointed within a few degrees to our line of sight. Such s mall viewing angles caus e ordinary 



29 



radio galaxies and radio loud quasars to appear as blazars (e.g.. lUrrv &: Padovanilll995l ). 



Optical variability can origi nate in relativistic jets or in the accretion disk (e.g.. lMarscher &: Gear 



19851 : iMangalam &: Wiitalll993l ). As the originating physical processes are different in each case. 



23 the light curves should look different. When FSRQs are in a high state, synchrotron emission from 

24 the relativistic jet overwhelms emission from the AD. Whe n FSRQs are in a low state, the AD can 

25 become visible, and be recognizable as the Big Blue Bump (jSun &: Malkanlll989l ). During these low 

26 and high states, which can differ by 2 to 5 magnitudes for the most active blazars (optically violent 

27 variables; OVVs), the dominant light source should determine the variability characteristics of the 
2B FSRQs, specifically, the power spectral densities (PSDs) and possible existence of quasi-periodic 

oscillations. 



30 The Kepler mission (iBorucki et al.ll2010l : iKoch et al.ll20ld ) is uniquely able to probe the inner- 

31 most regions of AGN and produce superior light curves through its ability to monitor uninterrupted 

32 by long gaps; it has only brief one-day gaps every month for data downlink. This capability allows 

33 us to study a broad range of variability time scales. The overwhelming majority of monitored 

34 blazars display signifi cant m i crova riability (variations of at least 0.03 mag) on timescales less than 

35 a day (see review bv iMillerl 1 19961 ): however, it is not known if this behavior is characteristic of 

36 all blazars because observers tend to select the most variable objects to monitor. We selected 

37 radio-loud AGN in the Kepler field of view for variability monitoring observations obtained over 8 
3B quarters in 2010-2012. Variability probably probes the characteristics of the accretion disks when 

39 the sources are in quiescent faint states and allows us to study the synchrotron jet emission when 

40 the sources are in highly active states. For both radio-loud and radio-quiet quasars, there is a good 

41 probability of observing significant variabil i ty (a few tenths of a m agnitude) over the course of two 

42 years (e.g. iPica et al.lll988l : iHawkinsI 120021: iMacLeod et al.ll2012l ). In a study of blazars observed 



Bauer et al.l (120091 ) found that 35% of blazars showed V > 0.4 



43 with the Palomar Quest survey, 

44 magnitudes of variation over 3.5 years. 

45 In addition, Kepler is capable of detecting microvariability over the course of several hours 

46 at the level of a few percent for targets brighter than about 17th magnitude. Ground-based op- 

47 tical studies have indicated that low frequency peaked blazars and core-dominated radio galax- 
4B ies with high polarizations frequently exhibit stronger microvariability, as expected if the jets 
49 point toward us and benefit from Doppler boosting. Ordinary radio-loud quasars and radio-quiet 



quiet quasars also show some microvariability, bu t it is both less free 



51 (e.g.. iGopal-Krishna et al.ll2003l : ICarini et al.l 120071 : iRamirez et al.ll2009l : iGoval et al.ll2012l ). Core 



uent and typically weaker 
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52 dominated radio-loud quasars, which are beheved to have jets pointing close to our line-of-sight, 

53 and therefore expected to show substantial fast variations, exhibit much more microy ariability than 

54 radio-quiet quasars if they also exhibit high optical polarizations (jGoval et al.ll2012l ). 



This paper is the first report on Kepler data for FSRQs. iMushotzky et alJ (120 111 ) reported 
Kepler monitoring of 4 Seyfert galaxies {z < 0.09), each for a duration of 2 to 4 quarters; we 
note that our targets have larger redshifts and are substantially fainter (by 2 to 3 magnitudes). 
According to their analysis, all four Seyferts showed some degree of variability over these periods 
and all exhi bited very stee p red n oise components to their power spectral densities (PSDs; —2.6 > 
a > —3.3). ICarini &: Rvld (j2012l ) provided a more in-depth analysis of the Kepler data for one of 
those Seyferts, II Zw 229.015, and were able to combine Kepler data across 3 quarters because they 
had sufficient ground-based measurements to normalize the fluxes on different detectors of Kepler's 
camera during different quarters. This allowed them to extend the PSD to lower frequencies where 
they detected a flattening below a frequency corres ponding to ~44 days; their best fit to the data 
in the regime measured by IMushotzky et al.l (120111 ) led them to a shallower PSD slope of around 
—2.8. Our results on FSRQs show significantly flatter PSDs than either group; we discuss this in 
§6. 

In this paper, we first briefiy review in §2 the physical origins of the variability signatures that 
disks and jets can be expected to exhibit, then describe our target sample selection in §3. In §4 
we explain the steps needed to reduce Kepler data on quasars, as it differs significantly from the 
standard analysis used for the primary Kepler science goal of exoplanet transit detection. We cover 
the post-processing analysis, including Palomar Observatory imaging and photometry, in §5. In §6 
we discuss our results, and finally summarize our conclusions in §7. 



2. Variability Signatures of Disks and Jets 
2.1. Emission from Accretion Disks 



Fluctuations emerging directly from the surfaces of ADs or from coronae above them can 
produce rest frame variations no faster than about t^ar — GMbh/c^ or ~ 8 minutes for Mbh = 
10^ Mq. The fastest real variations can provide a way of measuring the lower limits on masses 
of the SMBHs. Most X-ray and optical variability of the best studied cases, Seyfert galaxies and 
radio-quiet quasars, is "red noise", i.e., where the PSD arising from the Fourier transform of the 
light curve is characterized by P{f) oc with a < 0, below some break frequency, /;,; above 
fb the PSD is usually dominated by measurement errors with a white noise character (a = 0) or 
perhaps by white noise with an astrophysical origin. The PSD can also yield the size of the largest 
emitting region i f there is a flatter slope , often a ~ — 1 at low frequencies and ~ — 2 at intermediate 



85 frequencies (e.g.. iMarkowitz et al.ll2003l ) before perhaps turning to ~ at the highest frequencies. 



Variations yielding red-noise can be produced if energy is released over a wide range of time-scales 
characteristic of different orbital periods and turbulent transport in ADs (e.g.. iMangalam &: Wiita 
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19931 ). 



2.2. Previous indications of quasi-periodic variations 

Recently, detections of quasi-periodic oscillations (QPOs) in five AGN have be en made, mostly 



91 in X- rays. The best case is that of the Narrow Line Seyfert 1 galaxy RE J1034+396 ()Gierlinski et al 



20081 ). where a period of a little more than 1 hour was detected. The presence of optical QPOs 



93 has been suggested for the high declination BL Lac object, S5 0716+714, for some time (e.g.. 



Quirrenbach et al. 



1991 



A wavelet analysis of archival optical spectra of S5 0716+714 indicated 
the presence of five nights with QPOs present at > 0.99 probabil ity, with central periods ranging 
between ~ 25 and ~ 73 minutes (jCupta. Srivastava &: Wiitall2009l). New observ ations of this same 
blazar led to an even stronger indication of an ~15 min QPO (jRani et al.ll2010l ). 



9B The simplest QPO model assumes that a single hot-spot in the inner portion of the AD is 

99 responsible, and that radiation from this disk is directly detected, e.g., when an FSRQ is in a faint, 

100 quiescent state. Given a measured period, P (in seconds) one can estimate the SMBH mass M via 



M 



3.23 X 10^ P 



(1) 



Mq (r3/2 + a)(i + 2)' 

101 where r is the hot-spot distance in units of GM/c?, a is the BH spin-parameter and z is the 



102 redshift (iGupta. Srivastava &: Wiital l2009l ) . The shortest periods are obtained for the innermost 

103 stable circular orbit, and the strongest radiation emerges near there. For low redshift objects, a 5 

104 hour period, which would be very difficult to detect from ground-based observations, but rather easy 

105 for Kepler to find, corresponds to SMBH masses of 4.0 X 10^ and 2.5 x 10^ M©, for Schwarzschild 

106 and extreme Kerr BHs, respectively. In practical terms, we need to observe for many months to be 

107 sensitive to ADs surrounding billion-solar-mass black holes because any large hot spots in ADs are 
lOB probably both rare and short-lived, lasting no more than dozens of orbits (~ 2 weeks). 



109 2.3. Emission from Jets 

no For core-dominated flat spectrum radio quasars, where the emission is almost certainly dom- 

111 inated by jets within several degrees to our line of sight, the radio-X-ray variability can result 

112 from changes in a variety of physical parameters: the synchrotron-emitting population of parti- 

113 cles; the conditions under which they move; the orientation to line of sight; or the Doppler factor, 

114 5 = [7(1 — /3cos0)]~^, with /3 = V/c the velocity of the shock through the jet, and 9 the angle 

115 between the observe r's line of sight and the jet axis.. Major flares arise from new shocks pass- 

116 ing through the jets (IMarscher k, Geaiill985l '). but smaller fluctuations could arise from turbulence 

117 behind those shocks (IMarscher et al.ll2008l ) or "mini-jets" (jGiannios. Uzdenskv &: Begelmanl l2009l . 



2OIOI : iNalewaiko et al 



2011 



Alternatively, variations in the magnetic fields, changes in the ambi- 
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ent medium, and changes in the injected particle distribution can also cause variations in emitted 
light. The light travel time across the radiating region, modified by the Doppler factor, is the char- 
acteristic time scale. If these variable emissions arise from shocks in relativist ic jets, as is expected 
to be the case for FSRQs in high states, then the observed t^ar — Ai?(5c)~^ (jGopal-Krishna et al 



20031 ). Here AR is the physical size of the emitting region in its rest frame (roughly a jet diameter) 



124 and 5 the Doppler factor. In this case, variability timescales can constrain the jet's 6 (jJorstad et al 



20051 ). Constraints on the jet velocity and viewing angle can be particularly tight if other data, 



such as apparent superluminal motions of radio knots, where Vapp = ^ sin 0(1 — (3 cos 9) ^, can 
be detected through VLBI. Our 4 FSRQ targets are all VLBA calibrator s with compact rad io 
structures (e.g., [HB89] 1924+507=4C50.47, a member of the CJF survey in lBritzen etaD (12003)). 



A leading model for any jet-based quasi-periodic var iations is the int ersection of an outward 
propagating shock with a helical structure within the jet (iRani et al.ll2009l ). Even small deviations 
in the viewing direction result in significant changes in flux. The period of a QPO and its stability 
between events in a given source is a mode discriminator. Helical structures in a jet are relatively 
long-lived, so that QPOs associated with them should have similar periods during succe s sive events. 



134 But if QPOs involve turbulent cells behind the shock in the jet ([Marscher et al. 



2008 



Rani et al 



20091 ). short-lived dominant eddies are expected to give rise to different periods. 



2.4. Duty cycles and timescales of variable emission from relativistic jets 
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Dramatic variability is a defining characteristic of blazars; however, the duty cycles of activity 
are still not fully characterized, since the number of objects subject to sustained monitoring from 
the ground is not large. FSRQs vary on timescales of minutes to decades, with the most variable 
(e.g., 3C 279) rangin g two magnitudes from their a verage level and up to five magnitudes above their 
low quiescent level (iKartaltepe &: Balonekl 120071 ). Microvariability, at the level of 0.03 mag over 
the course of a night (typic a lly 6 hours) is quite common, with du ty cycles around 50% for blazars 
(|Gopal-Krishna et al.l l201ll : ICarini et al.l 120071 : iGoval et al.l |2012| . and references therein) . When 
the smallest variations detectable from the ground, at about 0.01 mag per night, are included, the 
blazar duty cycle rises to ~ 80%. The duty cycle for similar small variations from normal radio loud 
(steep spectr um) quasarSj for which jets are present, but not strongly Doppler boosted, is closer 



to 20% (e.g., ISagar et al. 



2004 



Stalin et al.ll2004l : iRamirez et al.ll2009l ). which might be expected. 



More surprisingly, radio-quiet QSOs show a similar duty cycle to radio-loud quasars, indicative of 
varia bility arising from accretion disks or perhaps fro m jets that are present only on nuclear scales 
(e.g.. iGopal-Krishna et al.ll2003l : iRamirez et al.ll2009l ). For one example, CTA 102 is a well s tudied, 
but otherwise typical, highly polarized bright {V ~ 15) blazar. lOsterman Meyer et al.l (120091 ) found 



variations as large as 0.1 mag in 15 minutes, with typical slopes of 0.01 mag/hr. 

Structure function analysis is used to characterize the variability timescale and infer the 
physical scales when the Doppler factor and orientation to our line of sight are also measured. 



Bachev et al. I (l2012l ) monitored the short term optical variability of 13 highly variable blazars (not 
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156 a complete sample) over the course of five years. They found surprisingly low variability: several 

157 objects varied slowly and smoothly at rates of up to ~ 0.1 mag per hour, but many displayed no 

158 short term variability. They quote only a ~ 2% chance of observing variability of more than 0.1 

159 mag per hour during their observations. 



160 3. Sample Selection 

161 We selected objects with strong compact cores by searching the NRAO VLBA Calibrator Lis10 

162 for radio sources brighter than 100 mJy on '^400 km baselines at 2.3 GHz in the Kepler field of 

163 view. We found four AGN that met these criteria. The properties of these targets are shown in 

164 Table [U for convenience, we refer to them throughout this paper as Objects A-D. Three objects 

165 are classified as quasars and have radio spectral indices a between 0.0 and 0.32 (where S oc 

166 the fourth object is a tentative Seyfert 1.5 galaxy which has a = 0.75 , typical of steep spectrum 

167 (extended) radio emission. We also observed the archetypal powerful radio galaxy, Cygnus A, with 
16B various custom apertures coordinated by Michael Fanelli at the Kepler Science Operations Center. 

169 Analysis of Cygnus A requires special treatment due to its extended optical structure and crowded 

170 foreground stars; we expect to present data and results in a subsequent paper. We observed all 

171 five targets in Kepler Guest Observer Cycles 2 and 3, corresponding to Quarters 6-13 (two years 

172 of almost continuous monitoring in 2010-2012, Kepler GO Programs GO20018 and GO30010), and 

173 we report these results and our analysis here. Further data on our targets (Cycle 4) is currently 

174 being acquired. 



175 4. Data Reduction and Analysis 

176 4.1. Overview 

177 The Kepler focal plane has 21 modules each with two 2200x1024 pixel CCDs for astronomical 

178 observing and 4 CCDs for fine guidance. The images are out of focus to improve sensitivity to 

179 periodic signals for extrasolar planet detection around bright stars. Thus the signal from even a 

180 point source is spread over several 2.98" pixels. The signals from pixels in an aperture whose siz e 



181 is a function of anticipated magnitude are sent down by the Kepler spacecraft (jBryson et al.ll2010l ). 

182 For faint targets such as our AGN (V magnitudes 17.8-18.6), the signal is several hundred electrons 

183 per second, a far different regime for instrumental response (and artifacts) than for the bright stars 

184 (y < 10) searched for extrasolar planets. We observed all of our targets in Long Cadence (LC) 

185 30-minute mode, yielding ~4500 samples each quarter per target. In addition, we observed our 

186 targets in Short Cadence (SC) 1-minute mode for a single quarter each year (~ 125,000 samples). 

187 Object B falls on dead CCD module #3 for one quarter of each year (Q8 and Q12 during Cycles 



http://www.vlba.nrao.edu/astro/calib/ 



-7- 



188 2 and 3), so we have only 6 quarters of data for it. In LC mode for Object D at V = 18.4 in 

189 the Kepler Input Catalog, the measured per-cadence noise level (Icr) was 0.9%. In SC mode, the 

190 per-cadence noise level at V = 18.4 was 3.4%. A journal of the observations is given in Table 2. 

191 The standard Kepler pipeline data have the day-to-week-scale drifts removed because the 

192 pipeline is optimized for extrasolar planet detection, which also removes real astrophysical bright- 

193 ness variations of the same timescales. Therefo re, we used pre-pipeline data, "SAP FLUX" , in our 



194 analysis. This was also the approach taken bv iMushotzkv et al 



(|201lh and ICarini k Rvld ([2012) 



195 in their analyses of Kepler data for brighter AGN. This choice means that some of the long term 

196 variations that are seen in the data are instrumental and not intrinsic to the source, but we can 

197 also be confident that we are not discarding any of the short term data trends that are intrinsic 

198 to the source. We removed the largest contribution of this effect on the PSDs we computed by 

199 performing end-matching (discussed below) on each quarter (for LC data) or each month (for SC 

200 data) . 



201 4.2. Instrumental Effects: Differential Velocity Aberration, Aperture Optimization, 

202 and Nonlinear Amplifier Electronics 

203 We note that there is no precise way to translate the observed signal in electrons/second to 

204 optical magnitude because this was not a mission design feature; Kepler was designed to measure 

205 very precise changes in relative brightness, not very precise absolute brightness. The light curves 

206 formed by Kepler long term data are affected at the few percent level by differential velocity 

207 aberration which originates in the motion of the spacecraft in its annual orbit around the sun. As 

208 the spacecraft moves, it keeps a constant field of view in the constellation Cygnus, and the angle 

209 formed by its velocity is constantly changing. This causes a target's PSF to move with respect to 

210 the detector array, and because a limited number of optimal aperture ("postage stamp") pixels is 

211 downloaded, each quarter's light curve would have a characteristic rising or falling slope or a concave 

212 or convex shape, superposed on other effects such as thermal drift. The finite number of pixels 

213 used for photometry is also responsible to the steps seen between quarters; after a spacecraft roll, 

214 an object falls on a different CCD and in a different location with respect to the pixel array. Five of 

215 the 16 CCDs that our four AGN land on each year are affected by an unstable amplifier electronics 



216 prob lem which causes a quasi-periodic signal to move across the CCDs (see iKolodziejczak et al 



2010 for technical details). This Moire effect, so-named because of its resemblance to a Moire 



218 pattern, affects faint sources more than bright sources because it is additive, not multiplicative. As 

219 of writing, no models had been developed to remove the effect from faint source signals. Affected 

220 CC Ds and modules are list ed in Table 13 of the Kepler Instrument Handbook; see also their Figure 

221 24 (jKepler ProiectI l2009al). and for mo dule number identification, see Figure 2-1 of the Kepler 

222 Archive Manual, (jKepler Proiectll2009bl ). These "Moire patterns" look like ripples with timescales 

223 of days, and are visible in data from some quarters, as listed in Table 2. The vast majority (80%) 

224 of our data are of excellent quality and are not affected by this problem. 
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225 4.3. Light Curves 

226 We show hght curves for ah four AGN in Figure 1. The time axis is given in Barycentric Julian 

227 Date (BJD), the standard Kepler time unit which is referenced to the Solar System barycenter; 

228 it differs from the more familiar (heliocentric) Julian Date by ±4 seconds (the recently discovered 

229 absolute timing error in the FITS headers (M. Still, Kepler Blog of 12 Dec 20lfl) has no effect 

230 on our studies). The repeated annual patterns of differential velocity aberration are the clearest 

231 features, producing long term smooth variations. 



232 5. Light Curve Analysis and Results 

233 5.1. Flares in CGRaBS J1918+4937 and MG4 J192325+4754 

234 We found a ~ 7% rise in the flux of CGRaBS J1918+4937 (object C) during QIO which lasted 

235 about five or six days around BJD=2455755 (Fig. 2). We confirmed that there was no plausible 

236 contamination source, e.g., a nearby star that might have flared and been incorporated into the 

237 pixels included in its photometry. This variation shares no characteristics with known artifacts, 
23B hence, we conclude that it is real. Object C also showed three consecutive fluctuations of ~ 3% 

239 amplitude during a span of about 25 days during Q9. Other variations in Object C during Q7 

240 and Qll might be real, but because there was some Moire effect present in those quarters, we 

241 cannot make that attribution. Object A (MG4 J192325+4754) showed several significant flares 

242 of 4 — 6% size over a period of 10 days in Q12. This quarter is one of the two for that AGN 

243 that should not have b een affected at all by Moire problems but this quarter was affected by solar 



244 coronal mass ejections (jKepler Proiectll2012l ): however, the interval in question was well away from 

245 those disruptions and also well away from the monthly gaps and glitches. Neither object B nor D 

246 showed any short-term variability of any type during any quarter. 



5.2. Search for Brief Flares 



No fast astrophysical flares (shorter than a few hours) were detected in any of our targets. 
Our search method was as follows. Two of us indep endently searched the S C and LC light curve 
data for flares. We used the Pyke task Key yixseries (jStill and Barclavll2012l ) to examine the data 
in each target's pixel; then we used topcat ( Taylor 20051 ) to find the flag values for anomalously 
high data points, identified in the Kepler Archive Handbook, Table 2-3. Both the two significant 
(defined at 4 to 6 consecutive high points each) flares turned out to be instrumental and not as- 
trophysical. The problems were flagged iii the F ITS flies as "coarse pointing" but were actually 



255 probably "argabrightening" (j Jenkins et al.ll20ld ) that was not recognized as such by the pipeline 



^http://keplergo. arc. nasa.gov/Blog.shtml 
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256 software. During Q12, strong solar coronal mass ejections caused spacecraft coarse pointing prob- 

257 lems that wer e manifested as man y anomalously high data points along with several anomalously 
25B low readings (jKepler Proiectll2012l ). 



5.3. Power Spectral Densities 

We constructed PSDs for all the Long Cadence data, then fitted slopes to the low frequency and 
high frequ ency portions of those PSDs. We computed PSDs with a Discrete Fourier Transform using 



262 Period04 (jLenz Sz Bregeij 120091 1 . both for the original SAP count rates and our "corrected" data, 

263 which underwent the following processing steps that fixed bad data points and would also allow 

264 us to use a fast Fourier transform on the light curves (see Fig. 2). First, monthly "glitches" were 

265 removed. These instrumental effects from the temporary changes in CCD sensitivities produced 

266 by the temperature changes resulting from reorientation of the spacecraft required to send back 

267 data each month often were at the level of 2-5% of the flux. We modeled them as exponential 

268 decays, with amplitudes and time constants as free parameters and then subtracted (or added) 

269 the corrections needed to produced smooth post-glitch curves. The typical decay time was 0.8 to 

270 1.5 days. Second, the gaps produced by the monthly data downloads, along with the occasional 

271 other gaps in the data (usually just a few dozen isolated points per quarter, occasionally two to 

272 five contiguous points) were filled in by fitting third order polynomials to the data on either side 

273 of the gaps and then adding noise derived from the local standard deviation. In the case of Q12 

274 data, which was affected by coronal mass ejections, these additional gaps lasted up to 4 days. 

275 Third, single point outliers more than 4(T from the local mean were removed and replaced with 

276 points generated from that local mean with a random dispersion based on that a value. Finally, in 

277 order to minimize errors in the PSD produced by the low-frequency instrumental drifts discussed 
27B above and to take into account t he fact that we do not have an infinit ely long data train, we "end- 



matched" each light curve (e.g.. lFougerdll985l : iMushotzky et al.ll201ll ). removing a linear trend so 



that the last several points are at the same level as the first several points in each quarter. 

The PSD slope is determined from equally weighted points in the PSD. The PSD is not 
uniformly populated with data; there is far more data at the high frequency end than the low 
frequency end. The slopes at the high frequency end are invariably consistent with 0, or white- 
noise, and our fits were done so as to force the high frequency PSD to a zero slope. Then the low-end 
PSD slope was found via a least squares procedure computed starting at the lowest frequency points 
until that slope started to fiatten dramatically, which was indicative of the location of the break 
frequency. The results are summarized in Table 3, which includes results for the original SAP data, 
end-matched SAP, corrected SAP and corrected-end-matched SAP. The slopes based on original 
SAP counts are usually between —1.5 and —2.0, but vary somewhat from quarter to quarter for 
the same source and from source to source, with Object D usually having a shallower slope than 
the other three. The corrected data yielded very similar slopes to those of the original data, 
though they were usually slightly steeper. It is very difficult to estimate the actual errors in the 
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293 PSDs (see lVaughn et al.ll2003l ). but the quarter-to-quarter differences between slopes of the same 

294 object are usuaUy characterized by standard deviations of 0.1-0.2, so the uncertainties involved in 

295 the computations should be no larger than this, since there is probably some degree of intrinsic 

296 variation between quarters. 

297 The slopes found using end-matched SAP data are almost always shallower than those using 
29B the pure SAP counts, which is expected since some of the long-term drift, which yields power at 

299 the lowest frequencies, is removed with end-matching. They are usually between —1.2 and —1.9. 

300 Note that if the intrinsic PSDs were actually steeper than —2.0, ignoring end-matching would tend 

301 to drive the estimated slopes toward around —2.0, while incorpo rating it provides much closer 



302 estimates of the actual slopes fed into simulations (jFougerd Il985l ). The fact that end-matching 

303 typically slightly flattens the slopes indicates that the intrinsic slopes were not steeper than —2.0. 

304 We show typical examples in Figure 3. 

305 We also fitted the Short Cadence data and give, in Table 4, the results for the original SAP 

306 data, the corrected SAP data and the corrected, end-matched SAP data; however, the last two 

307 values were always the same to at least two significant figures. For each quarter for each source, 
30B the break frequencies in the SC data turned out to be essentially the same as in the LC data, which 
309 reinforces our confidence in the fitting. 



310 5.4. Palomar Observatory Photometry 

311 Blazars can be "bluer when brighter" or "redder when brighter". Accordingly, we obtained 

312 optical and near-infrared photometry from CCD observations of our targets at the Palomar Obser- 

313 vatory 60 inch telescope in September, November and December 2010, during Kepler Q7. We used 

314 Sloan g'r'i' z' filters to evaluate the overall spectral shape within the very broad Kepler bandpass 

315 (Kepler's CCDs, which have no filters, are sensitive to visual wavelengths between from 4200 — 9000 

316 A). Data on the target AGNs and five to six nearby comparison stars in each field were reduced 

317 with the aperture photometr y tasks in IRAF. D ata were calibrated using the stellar magnitudes 



318 in the Kepler Input Catalog (iBrown et al.ll201ll ). Dereddening was applied using values from the 



319 NASA Extragalactic Database (NED). Results are given in Table 5. Objects A and B have nearly 

320 power-law spectra consistent with nonthermal emission, with a = 1.5 and 2.9, respectively for 

321 S oc z^~"). For comparison, blazar 3C 454.3, observed on 23 Nov 2010 in the course of another 

322 Palomar program, had a spectral index a = 2.0. In contrast. Object C has a nearly fiat spectrum 

323 (bluer than the other three o bjects) consistent w ith a blue continuum spectrum showing a single 



emission line, Mgll at 2798A, ( Healev et al.ll2008l : FITS spectrum provided by M. S. Shaw, personal 



325 communication). In Object C, therefore, the optical emission could have a significant contribution 

326 from thermal Big Blue Bump emission, but disentangling the relative contributions of nonthermal 

327 and thermal emission requires high quality spectrophotometry which is beyond the scope of this 
32B paper. Object D has a spectrum peaked at i' band, but the large uncertainty at z' band is sufficient 
329 to encompass a power law similar to those seen in Objects A and B. 
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5.4-1- Upper Limits on Variability 

None of the four AGN varied significantly during the Palomar observations in September 
through December 2010. At r' band, the brightness of the four targets was similar to those ob- 
served at the 1.2-m telescope at Mt. Hopkins in 2003-2008, as listed in the Kepler Input Catalog 



334 (jBrown et al.l 120111 ). Differences of ~ 0.1 — 0.4 magnitudes in three targets at g' , i' and z' were ob- 

335 served between the Kepler Input Catalog (KIC) measurements and the Palomar measurements; the 

336 reddest AGN (whose spectrum is most similar to the red stars used for calibration) showed no dis- 

337 crepancy. Three quasars are always fainter in blue bands and brighter in red bands than they were 

338 in the KIC. We attribute the Palomar-KIC differences to the effects of reddening and CCD-filter 

339 responses to the dissimilar spectra of quasars vs. stars. Observations at z' , where the discrepanc y 



is the largest, are most affected by site- variable atmospheric water vapor (jFukugita et al.lll996l ). 
moreover, different CCDs have different z' responses due to sensitivity variations at the longest 
wavelengths. 



5.4-2. Morphology 



Three of the four targets were unresolved in all four Sloan bands. Object B, the lowest redshift 
target, was resolved by ^'-band imaging into a central point source embedded in diffuse emission 
typical of host galaxies. With ground based seeing of l."5 and with the target at redshift 0.513, no 
structural details were expected or revealed; the fuzzy host galaxy emission in the i^'-band image 
was estimated at approximately 6 arcsec extent using ZR^F tasks. In the other Sloan bands, Object 
B was much more strongly dominated by the central nuclear source with comparatively weak diffuse 
emission detected in radial profile plots. It was tentatively clas sed as a Seyfert 1.5 ga laxy based on 
observations of a spectrum with a single emission line, Ha, bv lHenstock et al.l (j 19971 ). 



352 6. Discussion 

353 6.1. Flares and the Overall Shape of the Light Curves 

354 In 30 quarters of data, the strongest flare we saw was an excursion of ~ 7% amplitude in 

355 Object C during QIO, and lasting about 5 days. Object A also several significant flares of ~ 4 — 6% 

356 size over a period of 10 days in Q12. Object C showed three excursions of ~ 3% in Q9. Objects B 

357 and D showed no days-scale variability during any quarter. The flares could be caused by magnetic 

358 reconnection events or random fluctuations within the jets or by brightening in the accretion disk. 

359 No microvariability on timescales of minutes to hours was detected, possibly because the short- 

360 cadence noise levels of 2-3.4% were too high; ground-based observations show levels of 0.01-0.03 

361 magnitudes (~ 1 — 3%) are common in blazars, but rare in other classes of AGN, as reviewed in 

362 our Introduction. 
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We did not detect any astrophysical QPOs, which was unsurprising given the rarity of such 
events. Given the very long, consistent, Ught curves produced by Kepler, we greatly improved the 
chance of finding QPOs, if they were to be present on timescales of tens of minutes to weeks and 
sufficiently strong. We quantified our sensitivity to sine waves injected at a range of frequencies 
and amplitudes: in a blind study where one of us added such signals to real data another of us 
discovered that he could detect them in the PSDs at a level of ~ 3% of the original amplitude for 
18th magnitude targets, but could not do so for periodic signals at 1% amplitude. 

The com parative smo o thnes s of the Kepler light curves within each quarter is consistent with 
the finding by I Bauer et al.l (|2009l ) that only 35% of blazars showed V > 0.4 magnitudes of variation 
over 3.5 years, as observed with intermittent sampling during the Palomar Quest survey. Longer 
term optical variability of quasars has traditionally involved frequent monitoring of a modest num- 
ber of objects from individual observatories; those observations indicate that essentially all quasars 



375 do vary substantiall y over decades, and changes over a year or two are expected (e.g., iPica et al 



1988 



Hawkinal2002l ) . Averages of such monitoring c an yield statistical measures of quasar variabi L 



ity through the use of structure functions (SFs, e.g. Hughes et al. 1992 : Collier Sz PetersonlboOll ). 
Such monitoring, however, can be biased in favor of "interesting" objects, with those known to be 
variable being covered more frequently and some "uninteresting" objects can be dropped from the 
programs in favor of newly discovered variable objects. With the advent of massive sky surveys 
that revisit patches of the sky, such as the Sloan Digital Sky Survey (SDSS. lYork et al.ll200Cll ). a new 
approach has become possible: one can compute average variability estimates for tens of thousands 
of quasars divided into substantial bins by t heir redshifts, luminosities and other properties (e.g. 



Vanden Berk et al.ll2004l : Ide Vries et al.ll2005l ). Recently, structure function (SF) analyses of quasar 



variability using SDSS and the Palomar Sky Survey have led to estimates of quasar variability and 
a da n iped random wal k model that seems to explain its main properties (jMacLeod et al.l 120101 . 



20121 : iRuan et al.l l2012l ) . One conclusion of that work particularly relevant to our study is that 
radio-loud quasars are slightly more variable than radio-quiet ones. Specifically, the average SF^^f, 
a measure of the strength of variability, was arou nd 0.26 mag (in q') fo r the whole sample while the 
average for the radio- loud sample was 0.34 mag (jMacLeod et al.ll2010l ). where the redshifts in both 
categories were quite similar. Our densely sampled Kepler data on four AGN complement those 
programs with baselines of years that looked at very large numbers of objects much less frequently. 



393 6.2. Power Spectral Densities 

394 The results of our power spectral density (PSD) analyses are shown in Table 3 (long-cadence 

395 data) and Table 4 (short-cadence data). As discussed in §5.3, possible systematic errors in slopes 

396 are minimized by the corrections for (post-monthly-download) thermal drift and end-matching 

397 (last column of Table 3). The power-spectra for our quasars are dominated by red noise at lower 
39B frequencies and white noise at higher frequencies. We interpret the bulk of the red noise as intrinsic 
399 to the quasars; the white noise is instrumental. We find mean slopes between —1.2 and —1.8 for the 
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red-noise portion of the PSDs for our four AGN. These sfopes varied from quasar to quasar and, 
to a somewhat lesser extent, from quarter to quarter for the same object. Object D displayed the 
flattest slopes on the whole, with those for Object B were somewhat steeper but distinctly shallower 
than Objects A and C. Object B, tentatively classed as a Seyfert 1.5 galaxy, is the object in which 
our Palomar imaging detected the faint host galaxy surrounding the bright central point source, 
hence, the Kepler-detected optical emission had a small extra contribution from the host galaxy 
which could damp down any variation and flatten the resulting PSD. Among our four objects, there 
was no trend for PSD slopes to vary with redshift. The observed PSD s lopes are consistent with 



those of turbulence in sho cks in a relativistic jet or in an accretion disk (jMarscher &: Travis! 1 1991 



Mangalam fc Wiitalll993l ). Models of relativistic tur bulence which may apply to our quasars are 
being developed in the context of gamma ray bursts (|Zrake and MacFadvenll2013l ) and references 
therein. 



Kepler measurements of the brighter Seyfert galaxies analyzed bv lMushotzkv et al.l (j201ll ) have 
substantially steeper slopes (between —2.6 and —3.3) for the PSDs of their 4 Seyferts, and they also 
employed end-matching. Their objects were all substantially brighter than our targets (by factors of 
about 1.5 to 15) and so they could easily detect smaller fractional variations of 0.1-1%. Given that 
our objects are at much greater distances, and therefore substantially more powerful, despite their 
lower fluxes, there might be some physical reason for the discrepancy. Their results of very steep 
PSD slopes were surprising, because long-term X-ray m easurements of other Sey f erts yielded PSD 



slopes hat are almost always shallower thari —2.0 (e.g. lEdelson &: Nandr all 19991 : iMarkowitz et al 



Lachowicz et al. ! 


2009 


). 


Mushotzkv et al.l 


(2011 



20031 : lMarkowitzll2009l : 

ZW 229.015, for which iMushotzky et al.l (|2011l ) found slopes between —2.96 and —3.31, a reanalysis 
by a different method, coupled wi th ground-based observations yielded a best-fit shallower slope 
around —2.83 (ICarini &: Rylell2012l ). Our independent analysis of data for that object following the 
correction an d end-matching procedur es used in this paper produced slopes consistent with —2.2. 
We note that iMushotzkv et al.l (l201ll) remark that "PSD analyses are notoriously susceptible to 
analytical systematics (see, e.g. JVaughn et al.l (j2003l ))". Users of PSDs are on firmer ground when 
comparing "apples to apples", that is, PSDs for different objects are internally consistent when 
computed by each group's analysis techniques, but PSDs generated by different groups' analysis 
techniques are much more difficult to compare. 

The break frequencies at which the red noise and white noise contributions are equal were 
between ~ 10~'^'^Hz and 10~^'^Hz, corresponding to 0.2 to 0.6 days. Increasing the instrumental 
noise raises the white noise "floor" in the PSD diagrams, shifting the intersection of red noise and 
white noise to longer times cales (to the left i r i our figures). Kepler measurements of the brighter 
Seyfert galaxies analyzed bv lMushotzkv et al.l (j201ll ). which have smaller fractional noise, also have 
break frequencies of about 0.25 days. 
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436 6.3. Future Work 

437 We are currently acquiring additional long cadence data on each target. In future work, we 

438 will apply principal components analysis and cotrending basis vectors to the long cadence data to 

439 separate out common instrumental long term trends (such as differential velocity aberration and 

440 thermal drift, not Moire effects) in the data. These techniques may allow us to patch together 

441 individual quarters, forming light curves two or more years long for three targets and several three- 

442 quarters light curves for the remaining target that falls on a dead CCD during one quarter each 

443 year. Our goal is to detect lower b reak frequencie s in th e PSDs to find the largest physical scale 



444 as was done for II Zw 229.015 by ICarini &: Ryld (120121 ) who used ground-based data to bridge 



445 quarterly amplitude jumps in their Kepler light curves. 



446 7. Summary 

447 We have observed three flat spectrum radio quasars and one radio-loud AGN (tentatively 

448 classed as a Seyfert 1.5 galaxy) for two years with Kepler. We find power spectral densities in the 

449 light curves that vary from source to source and quarter to quarter for the same source, ranging 

450 from a = —1.2 to —1.8. We find several clear, isolated flares over a few percent amplitude which are 

451 of astrophysical origin. We do not detect any of the expectedly rare quasi-periodic oscillations as 

452 could be generated by accretion disks or helical jet features. The power spectral densities measured 

453 agree with models for the observed variability originating in turbulence behind a shock in the jets 

454 or in the accretion disks. These observations provide beautiful sets of high quality data. Future 

455 work on combining data across quarterly boundaries and extending the observations for the life of 

456 the Kepler mission could provide us with estimates of the largest relevant scales of the physical 

457 processes producing the variability. With the further development of shocked jet models including 

458 turbulence (e.g., Marscher, in preparation), these types of continuous light curves may distinguish 

459 between jet turbulence and accretion disk fluctuations as the dominant source of quasar optical 

460 variability. 
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Table 1. Kepler AGN Monitoring Target List 



Object 


Name 


Kepler Input 


Right 


Declination 


Kepler Input 


Redshift 


Radio 


Designation 




Catalog 


Ascension 




Catalog 




Spectral 






Number 






Magnitude 




Index" 








(hh:mm:ss.s) 


(dd:mm:ss) 








A 


MG4 J192325+4754 


10663134 


19:23:27.24 


47:54:17.0 


18.6 


1.520 


0.32 


B 


MG4 J190945+4833 


11021406 


19:09:46.51 


48:34:31.9 


18.0 


0.513 


0.75 


C 


CGRaBS J1918+4937'' 


11606854 


19:18:45.62 


49:37:55.1 


17.8 


0.926 


0.00 


D 


[HB89] 1924+507 


12208602 


19:26:06.31 


50:52:57.1 


18.4 


1.098 


0.19 



''Radio spectral index obtained from VLB A Calibrator website, defined between 2.3 and 8.3 GHz or 2.3 and 8.6 GHz with S oc u 

'^Kepler Input Catalog incorrectly indicates that this target is a star with contamination 0.73, but we have verified it is an isolated 
quasar. 



Table 2. Journal of Kepler Observations of AGN 



Object Long Cadence" Short Cadences'^ 

(Kepler quarters)'' (Kepler months) 

A Q6*, 7**, 8, 9**, 10*, 11**, 12, 13** Q6.1, 6.2, 6.3, 11.1, 11.2, 11.3 

B Q6, 7, 9, 10, 11, 13 Q6.1, 6.2, 6.3, 11.1, 11.2, 11.3 

C Q6, 7*, 8, 9, 10, 11*, 12, 13 Q6.1, 6.2, 6.3, 11.1, 11.2, 11.3 

D Q6*, 7, 8, 9, 10*, 11, 12, 13 Q6.1, 6.2, 6.3, 10.1, 10.2, 10.3 



"Targets are observed on a given Kepler detector for one quarter, using 4 detectors during a year. Start-stop dates 
for quarters are given at: http://keplergo.arc.nasa.gov/ArchiveSchedule.shtml For our Q6-13 data, each quarter 
returned approx. 4100-4700 30-minute cadences. Object B falls on dead CCD #3 during Q8, Q12, etc. 

''Some data are affected by a Moire fringe-like time-dependent instrumental fluctuation in amplitude that is very 
hard to calibrate. * denotes likely fringing; ** denotes fringing obvious by inspection of the raw data. 



'^For our Q6, 10 and 11 data, each month returned approx. 40,000 1-minute cadences. 
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Table 3. AGN variability power spectral densities (PSDs) measured with long-cadence Kepler 

data 



SAP (raw) data Corrected SAP data 





Kepler 




Moire 


oiope 




olope 


Oiope 


D 6sign 3,t ion 


ID 




level 


(original) 


\GH Q- m dX CilGCl j 


(original) 


(end-matched) 


A 


10663134 


6 


Medium 


-1.9 


-1.7 


-2.0 


-1.6 


A 


10663134 


7 


High 


-2.1 


-2.1 


-2.3 


-1.3 


A 


10663134 


8 


None reported 


-2.0 


-1.9 


-2.1 


-2.0 


A 


10663134 


9 


High 


-1.7 


-1.3 


-1.6 


-1.4 


A 


10663134 


10 


Medium 


-1.5 


-1.6 


-1.6 


-1.4 


A 


10663134 


11 


High 


-2.0 


-1.7 


-2.1 


-2.0 


A 


10663134 


12 


None reported 


-1.8 


-1.9 


-2.0 


-1.8 


A 


10663134 


13 


High 


-1.9 


-1.7 


-1.8 


-1.8 


Mean'^ 








-1.8 


-1.8 


-1.9 


-1.7 


SD 








0.2 


0.2 


0.2 


0.3 


B 


11021406 


6 


None reported 


-1.9 


-1.6 


-2.0 


-1.4 


B 


11021406 


7 


None reported 


-1.8 


-1.3 


-1.9 


-1.2 


B 


11021406 


8 


None reported 










B 


11021406 


9 


None reported 


-1.5 


-1.7 


-1.8 


-1.8 


B 


11021406 


10 


None reported 


-1.8 


-1.0 


-1.5 


-0.9 


B 


11021406 


11 


None reported 


-1.6 


-1.7 


-1.8 


-1.3 


B 


11021406 


12 


None reported 










B 


11021406 


13 


None reported 


-1.8 


-1.6 


-1.9 


-1.7 


Mean 








-1.7 


-1.5 


-1.8 


-1.4 


SD 








0.2 


0.3 


0.2 


0.3 


C 


11606854 


6 


None reported 


-1.9 


-1.4 


-2.0 


-1.6 


c 


11606854 


7 


Medium 


-1.5 


-1.9 


-1.6 


-2.0 


C 


11606854 


8 


None reported 


-1.8 


-1.2 


-2.0 


-1.6 


C 


11606854 


9 


None reported 


-1.9 


-1.9 


-1.9 


-2.0 


c 


11606854 


10 


None reported 


-1.9 


-2.4 


-2.1 


-2.1 


c 


11606854 


11 


Medium 


-1.8 


-2.0 


-2.0 


-2.0 


c 


11606854 


12 


None reported 


-1.8 


-1.7 


-1.9 


-1.7 


c 


11606854 


13 


None reported 


-1.8 


-2.0 


-2.0 


-2.0 


Mean 








-1.9 


-1.8 


-2.0 


-1.8 


SD 








0.1 


0.4 


0.1 


0.2 


D 


12208602 


6 


Medium 


-1.8 


-1.5 


-1.9 


-1.3 


D 


12208602 


7 


None reported 


-1.4 


-1.2 


-1.5 


-1.0 


D 


12208602 


8 


None reported 


-1.0 


-0.8 


-1.2 


-0.8 


D 


12208602 


9 


None reported 


-1.9 


-1.6 


-2.0 


-1.6 


D 


12208602 


10 


Medium 


-2.0 


-1.4 


-2.0 


-1.4 


D 


12208602 


11 


None reported 


-1.9 


-1.3 


-1.8 


-1.4 


D 


12208602 


12 


None reported 


-1.5 


-1.2 


-1.8 


-0.9 


D 


12208602 


13 


None reported 


-1.1 


-1.2 


-0.7 


-1.3 


Mean 








-1.5 


-1.2 


-1.5 


-1.2 


SD 








0.4 


0.2 


0.5 


0.3 



^Means and standard deviations computed using only quarters with no Moire reported, except for A, for which Q6 and QIO 
were also included. 
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Table 4. AGN variability power spectral densities (PSDs) measured with short-cadence Kepler 

data 











SAP (raw) data 


Corrected SAP data 


Object 


Kcplcr 


Quarter 


Moire 


Slope 


Slope 


Designation 


ID 


-Month 


level 


(original) 


(original & end-matched) 


A 


10663134 


6-1 


Medium 


-2.0 


-2.0 


A 


10663134 


6-2 


Medium 


-2.1 


-2.0 


A 


10663134 


6-3 


Medium 


-1.7 


-1.9 


A 


10663134 


11-1 


High 


-1.6 


-2.0 


A 


10663134 


11-2 


High 


-2.0 


-2.5 


A 


10663134 


11-3 


High 


-1.6 


-1.6 


Mean'' 








-1.9 


-1.9 


SD 








0.2 


0.1 


B 


11021406 


6-1 


None reported 


-2.0 


-1.9 


B 


11021406 


6-2 


None reported 


-1.7 


-1.7 


B 


11021406 


6-3 


None reported 


-1.3 


-1.7 


r> 




11 1 

11 — i 


None reported 


-i.o 


1 O 

-L.Z 


B 


11021406 


11-2 


None reported 


-1.3 


-1.3 


B 


11021406 


11-3 


None reported 


-1.8 


-1.9 


Mean 








-1.6 


-1.6 


SD 








0.3 


0.3 


C 


11606854 


6-1 


None reported 


-1.8 


-2.0 


C 


11606854 


6-2 


None reported 


-1.8 


-1.8 


C 


11606854 


6-3 


None reported 


-1.9 


-1.8 


C 


11606854 


11-1 


Medium 


-1.8 


-2.0 


C 


11606854 


11-2 


Medium 


-2.0 


-2.0 


C 


11606854 


11-3 


Medium 


-1.9 


-1.8 


Mean 








-1.9 


-1.9 


SD 








0.1 


0.1 


D 


12208602 


6-1 


Medium 


-1.9 


-1.7 


D 


12208602 


6-2 


Medium 


-1.7 


-1.8 


D 


12208602 


6-3 


Medium 


-1.8 


-1.8 


D 


12208602 


10-1 


Medium 


-1.8 


-1.9 


D 


12208602 


10-2 


Medium 


-1.4 


-1.7 


D 


12208602 


10-3 


Medium 


-1.8 


-1.8 


Mean 








-1.7 


-1.8 


SD 








0.2 


0.1 



*Mean and standard deviations computed using only quarters with no or medium Moire reported. 
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Table 5. Palomax Observatory Photometry 



Object 


Filter 


Wavelength 


Frequency 


KIC 


Palomar 


Error 


Reddening 


De-rcddcncd 


Dc-rcddcncd 


Flux Density 




Band 














Palomar 


Flux Density 


Error 






(jiui) 


(ilz) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(!"■)>-) 


(m.)y) 




g' 


0.52 


5.77E+14 


18.597 


19.15 


0.02 


0.345 


18.80 


0.11 


0.02 


A 


r' 


0.67 


4.48E+14 


18.591 


18.62 


0.02 


0.244 


18.37 


0.20 


0.02 


A 


i' 


0.79 


3.80E+14 


18.652 


18.46 


0.02 


0.189 


18.27 


0.23 


0.02 


A 


z' 


0.91 


3.30E+14 


... e 


18.34 


0.11 


0.161 


18.18 


0.26 


0.10 


Bb 


g' 


0.52 


5.77E+14 


18.539 


18.88 


0.08 


0.229 


18.65 


0.13 


0.07 


B 


r' 


0.67 


4.48E+14 


17.973 


17.99 


0.05 


0.168 


17.82 


0.34 


0.05 


B 


i' 


0.79 


.3.80E+14 


17.790 


17.64 


0.06 


0.13 


17..51 


0.17 


0.06 


I! 


z' 


0.1)1 


.i.aOK+ll 




17. 2o 


0.1:5 


0.111 


17.11 


().()7 


0.12 


c<= 


g' 


0.52 


5.77E+14 


17.833 


17.91 


0.03 


0.288 


17.62 


0.33 


0.07 


c 


r' 


0.67 


4.48E+14 


17.723 


17.77 


0.01 


0.199 


17.57 


0.34 


0.05 


c 


i' 


0.79 


3.80E+14 


17.784 


17.93 


0.03 


0.148 


17.78 


0.28 


0.04 


c 


z' 


0.91 


3.30E+14 


17.617 


18.00 


0.26 


0.11 


17.89 


0.25 


0.25 


Dd 


g' 


0.52 


5.77E+14 


18.498 


18.84 


0.05 


0.383 


18.46 


0.15 


0.05 


D 


r' 


0.67 


4.48E+14 


18.394 


18.23 


0.05 


0.282 


17.95 


0.30 


0.04 


D 


i' 


0.79 


3.80E+14 


18.408 


18.15 


0.04 


0.215 


17.94 


0.32 


0.04 


D 


z' 


0.91 


3.30E+14 


... 6 


18.35 


0.29 


0.185 


18.17 


0.26 


0.27 



^Observed on 2010-Sep-22, 2010-Dec-05. 

''Observed on 2010-Nov-13, 2010-Dec-05, 2010-Dec-09. 

'^Observed on 2011-Nov-12, 2011-Dec-02. 

■^Observed on 2011-Nov-12, 2011-Dec-03. 

"No z' data were available in the KIC for these objects. 
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Fig. 1. — Kepler light curves for objects A-D (starting at top) covering Quarters 6-13 in 2010-2012. 
The vertical axes have units of electrons per second, the horizontal axes have units of days where 
BJD is the Barycentric Julian Date. See §4.2 for an explanation of the steps between quarters. 
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Fig. 2. — Kepler light curve for Object C in Quarter 10. The rise in the amplitude around BJD 
2455755 lasting five to six days is clearly visible, as are the two thermal excursions after the monthly 
data downloads. The vertical axes have units of electrons per second, the horizontal axes have units 
of days where BJD is the Barycentric Julian Date. 
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B Quarter 6 LC RAW Corrected 



B Quarter 7 LC RAW Corrected 




560 580 600 

Baiyccntric JD (Days -2455000) 



660 680 700 

Barycentric JD {Days -2455000) 



B Quarter 9 LC RAW Corrected 



i Quarter 10 LC RAW Corrected 



435 

430 
425 
420 
I 415 
410 
405 




840 860 880 900 

Barycentric JD (Days -2455000) 




940 960 
Barycentric JD (Days -2455000) 



450 1 

1000 



B Quarter 1 1 LC RAW Corrected 



i Quarter 13 LC RAW Corrected 



1020 1040 1060 1080 

Barycentric JD (Days -2455000) 



g 410 

< 400 



1200 1220 1240 1260 

Barycentric JD (Days -2455000) 



Fig. 3. — Kepler light curves for Object B for 6 quarters (no data were taken for Q8 and Q12 due 
to an inoperative detector). Several operations were performed on the "raw", or SAP, data prior 
to computing PSDs (see §5). Color coding represents: unmodified data (light blue); removal of 
monthly thermal glitch (orange); monthly data gap and other missing data (green); data entries 
with a time value but no flux value (purple); and outliers (red). The dominant variation is not 
intrinsic variability, but residual 1-year pe riod differential velocity aberration and thermal drift 
instrumental effects (jKinemuchi et al.ll2012l ). 
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Fig. 4. — Examples of flux variability power spectral densities (PSDs) for Objects A, B, C, and 
D. Each shows data from a single quarter, with fluxes corrected for known calibration issues as 
described in the text. A two-parameter fit is shown for each: (i) a white-noise level that dominates 
at high frequencies; and (ii) a power law slope that dominates at frequencies below ~ 10~^ Hz. 



